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Abstract 



There is discussion if traffic displays multiple phases (e.g. laminar, januned, synchronized) or not. This paper presents 
evidence that a stochastic car following model, by changing one of its parameters, can be moved from showing two 
phases (laminar and jammed) to showing only one phase. Models with two phases show three states: two being 
homogeneous states corresponding to each phase, and a third state which consists of a mix between the two phases 

(phase coexistence). 

Although the gas-liquid analogy to traffic models has been widely discussed, no traffic-related model ever displayed 
a completely understood stochastic version of that transition. Having a stochastic model is however important to 
understand the potentially probabilistic nature of the transition. Most importantly, if indeed 2-phase models describe 
certain aspects correctly, then this leads to predictions for breakdown probabilities. Alternatively, if 1 -phase models 
describe these aspects better, then there is no breakdown. Interestingly, such 1-phase models can still allow for jam 
formation on small scales, which may give the impression of having a 2-phase dynamics. 
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1 INTRODUCTION 

Both from an operations and from a design perspective, the capacity of a road is an important quantity. Clearly, if 
demand exceeds capacity, queues will form, which represent a cost to the driver and thus to the economic system. In 
addition, such queues may impact other parts of the system, for example by spilling back into links used by drivers 
who are on a path that is not overloaded. 

For a variety of reasons, however, capacity is not a deterministic fixed quantity. It is possible that one day a 
queue forms and the next day not, and this may even happen in spite of demand being larger on the second day. In 
consequence, any definition of capacity needs to take its stochastic nature into account. 

For example, one could measure flow in 15-min intervals, say from 6am to 6: 15am, from 6: 15am to 6:30am, 
etc. One could then take the daily maximum of these values, and average the result over many typical workdays. 

As an alternative, one could measure flow as a function of density irrespective of any other state variables. 
One could then obtain the average flow for each density interval, and the maximum of these flow-values would repre- 
sent capacity. 

All of these measurements have the property that they result in an expected value, i.e. in a number that, for 
a given day, can be exceeded or not be reached. In consquence, it is useful to develop models of traffic which reflect 
the stochastic nature of traffic. Clearly, the stochasticity can come at many different levels: demand can vary; road 
conditions can vary; driving behavior can vary; etc. These different contributions to stochasticity will have different 
influences, which need to be debated. In this paper, we want to concentrate on road capacity. We understand that there 
is active research to eventually include aspects of stochastic transitions into the Highway Capacity Manual (0). 

The starting point for our work are simple single-lane car following models. These models are typically either 
of the type v{t + t^) = f{g, Av, ...) or of the type a{t + Tq) = h{g, Av, ...), where v{t) is the velocity of a car at 
time t, Av is the velocity difference to the car ahead, and a is the acceleration, g is the gap to the car ahead, where 
g = Ax — I, with Ax the front-buffer-to-front-buffer distance, and I is the space the car occupies in a jam. These 
models can for example be found in 

v{t+T^) = V{g{S)) with V{g) =vf -Vf exp{-Xg/vf) , 

ind) 

a{t) = a ■ (v{g{t)) - v{t)^ , with V{g) = Vf ■ (tanh(3 + £) - tanh(^)) , 

or in 

^ Ht + Ta)y ^ 

Additional parameters here are Vf (the free speed). A, I, and m. 

When these models are implemented on a computer, they need to be discretized in time, and one has to 
concern oneself with the size of the integration time step. At. A typical discretization is 

a{t) = given by the model 

v{t) = v{t - At) + Ata{t) , 

and 

x{t + At) ^ x{t) + Atv{t) . 

These discretizations are meant to approach the original coupled differential equations for At 0, and there is a 
whole body of literature available for this. Once time delays (via t > 0) are introduced into such equations, numerical 
treatment becomes considerably more difficult, and much less is known about efficient numerical integration. 

Given that state of affairs, it makes sense to look for computational models which are not based on the limit 
At 0, but which generate useful results also for relatively large time steps of, say, one second. The model that 
we will use in this paper has been introduced by Krauss (^; it is a variant of a model used by Gipps (H). The Krauss 
model has been shown to be free of collisions (i.e. that g < never occurs). 

In addition to being crash free at large time steps, the Krauss model is also stochastic. The important param- 
eter for our study is a noise amplitude e, which we wiU vary from 0.5 to 2. For e < or e > 2 the model leaves the 
range of where it is plausible for traffic. 

Our main results are the following: 
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• For medium e, there are three states of traffic, which we will call laminar, mixed, and jammed. The state 
depends on the density. Laminar, occuring at low density, means that there are no stopped cars; jammed, 
occuring at high density, means that no car is driving at full speed; and mixed, occuring at intermediate 
density, means that the system is a mix of laminar and jammed traffic. In the mixed state, traffic is strongly 
inhomogeneous; the mixed state is often called the coexistence state. 

It is important to note that there are three states (laminar, jammed, mixed) but only two phases (laminar, 
jammed). The phases refer to homogeneous sections of the system; the state refers to the system as a whole. 

• For large e, there is only one phase of traffic and therefore only one state. When going from low to high 
density, cars move closer and closer together, but traffic remains homogeneous at all times. 

• At some e in between, there is a transition from the 1 -phase to the 2-phase regime. 

• In the Kxauss model, changes of e also change the average acceleration. This is an unfortunate coincidence, 
and we believe that our general results regarding the number of phases are not related to this effect. 

• Deterministic models, formulated either as car following models or as fluid-dynamical models, can display 
1 -phase or 2-phase behavior They can however not display stochastic transitions between the phases. 

The results are important for model building as well as for understanding field measurements. In a 2-phase 
model, theory predicts that there can be a hysteretic transition from the laminar to the mixed state without a change 
in density. This means that, at a given density, traffic can operate in the laminar flow state for long times, until it will 
eventually "break down" and switch to the mixed state. In a 1 -phase model, this is impossible, and there is only one 
phase for any given density. 

A direct consequence of this is that, if traffic follows a 1 -phase model, any initial jam will "smear out" 
and thus eventually go away, even with unchanged traffic conditions. Conversely, in a 2-phase model, jams "pull 
themselves together" and thus stabilize themselves. 

This paper starts with Sec. || which describes the general idea of a gas-liquid transition. Sec. [5] describes the 
general simulation setup including the car following model that is used, discusses space-time plots of the resulting 
dynamics, and investigates transients vs. the steady state. Sec. ^ then establishes how a mixed state can be nuemrically 
detected for a given model. Sec. ^ discusses how these results relate to deterministic model; the paper is concluded by 
a discussion and a summary. 



2 PHASES IN TRAFFIC 

At a first glance, our results may seem largely irrelevant to traffic operations. As already discussed above, we contend 
that it is not. There are thoughts about extending the Highway Capacity Manual so that it includes the concept of 
stochastic traffic breakdown (^. This could for example mean that, for certain flow levels, one would include a curve 
describing the probability that traffic flow has not broken down as a function of time. In order for such a description to 
make sense, the existence of spontaneous flow break-down needs to be established. Some researchers however have 
doubts that flow breakdown truly happens ^ |^). This paper will contribute to the theoretical understanding of the 
issue, by showing how one model can move from displaying breakdown to not displaying breakdown. 

In fact, the analogy between a gas-liquid transition and the laminar-jammed transition of traffic was pointed 
out many times (e.g. (^ ^). The description of traffic in the well-known 2-fluid-model ([To^ im plicitely assumes 
the existence of two phases; and all simulation models which use spatial queues (e.g. ( [77[ |/2[ [73|)) will display two 
phases because of the definition of the dynamics. The two phases in models with queues are however much easier to 
understand than the phases in more realistic models. 



In a gas-liquid transition, one observes the following (see also Fig. 1(a) left): 



In the gas phase at low densities, particles are spread out throughout the system. Distances between particles 
vary, but the probability of having two particles close to each other is very small. 

In the liquid phase at high densities, particles are close to each other There is no crystalline structure as in 
solids, but the density is similar and in some cases (e.g. in water) even high in the liquid than in the solid 
phase. Because of the fact that the particles are so close to each other, it is difficult to compress the fluid any 
further. 
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In between, there is the so-called coexistence state, where gas and liquid coexist. In typical experiments 
in gravity, the liquid will be at the bottom and the gas will be above it. Without gravity, as well as for 
example within clouds, droplets form within the gas and remain interdispersed. In clouds, small droplets will 
eventually merge together into bigger droplets, which will fall out of the cloud as rain. Without gravity, the 
droplets will just merge but never fall out. The final state of the system is having one big droplet of liquid, 
surrounded by gas. 



The kinetics of the droplet formation (e.g. (|74[)) is ruled by a balance between surface tension and vapor 
pressure. Since surface tension pulls the droplet together, it has a tendency to push particles out of the droplet. Vapor 
pressure is the balancing force - it pushes particles into the droplet. 

Surface tension depends on the droplet radius - the smaller the droplet, the larger the surface tension. The 
result is that, when coming from small densities, there is a regime where large droplets would already be stable, but 
small droplets are not. In equilibrium, this density would be in the coexistence state, but coming from low density, the 
homogeneous gaseous phase can survive for some time. This way, meta-stability can be explained in the gas-liquid 
transition, and by extension also in traffic models. 

When compressing the system beyond this point, at density pi, the following can be observed: Pressure will 
not go up any further. Instead, more droplets will form, and the density outside the droplets will remain constant. 
Pressure will not go up until all the space is used up by droplets, i.e. all particles are moved into the liquid phase. Only 
at this point, at density p2, will pressure go up again when compressing the system. 

The state between pi and p2 is called coexistence state. In it, gas and liquid coexist, and changes on the 
global density are absorbed by the fraction of space and particles that the liquid occupies. The interface between 
droplets and gas remains the same no matter what the global density is. 

In addition to the coexistence, when keeping density constant one will observe coagulation as already men- 
tioned earlier: The many original small droplets will coagulate to fewer and fewer large droplets; in the limit there will 
only be one large droplet. 

A direct consequence of the above behavior is hysteresis: 

• When coming from low densities, it is possible to go beyond pi and still remain in the gaseous phase. Only 
after some waiting time will, by a fluctuation, some particles get close enough to each other to start the 
formation of a droplet. 

• When coming from densities above pi, it is possible that the droplet survives for some time even below pi. 

When increasing temperature with the above model, the 2-phase structure will eventually go away. This 
happens via pi and p2 approaching each other and eventually meeting. That is, depending on the temperature, a 
fluid system will eith er di splay transitions from gas to coexistence and from coexistence to liquid, or there will be no 



transition at all (Fig. |l(a)| right). 

We will now move on to describe the supporting evidence for our claims. As is typical in computational 
statistical physics, our evidence is based on computer simulations. It is however backed up by generic knowledge 
about phase transitions as they are well understood in physics. 

3 THE SIMULATIONS 
3.1 Krauss Model 

The velocity update of the Krauss model ( [73| [76[ ) reads as follows: 

, gft) - i{t)T 

Wsafo = V{t) + _,, , (1) 

v/b + T 

Wdos = min{i;(<) + a Ai,i;safo,Wmax} (2) 
t'(t + At) = max{0, Wdcs ~ e a J?} ■ (3) 

V is the speed of the car in front, v — (w + w)/2 is the average velocity of the two cars involved, w^ax is the maximum 
allowed velocity, a is the maximum acceleration of the vehicles, h their maximum deceleration for e = 0, e is the noise 
amplitude, and 77 is a random number in [0, 1]. The meaning of the terms is as follows: 
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• Eq. (|l]): Calculation of a "safe" velocity. This is the maximum velocity that the follower can drive when she 
wants to be sure to avoid a crash. The main assumption is that the car ahead will never decelerate faster than 
b, and that the car of the follower can also decelerate with up to b. 

• Eq. (||): The desired velocity is the minimum of: (a) current velocity plus acceleration, (b) safe velocity, (c) 
maximum velocity (e.g. speed limit). 

• Eq. (|^: Some randomness is added to the desired velocity. 

After the velocities of all vehicles are updated, all vehicles are moved. 

The deterministic limit e = of the Krauss model has been proven to be free of crashes for numerical time 
steps At smaller than or equal to the reaction time, t. We will use At = r = 1 as has conventionally been used for 
the Krauss model. Even for e > 0, as was used for our studies, we never observed vehicles getting closer than their 
minimum distance. We further use a = 0.2, b = 0.6, v,nax ~ 3 for all simulations. 

The model is free of units; this is a property that it has inherited from the cell-based cellular automata models. 
A reasonable calibration is: time steps correspond to seconds and cells correspond to 7.5 meters. The reaction time 
then is assumed to be 1 sec, and u^ax — 3 corresponds to 22.5 m/s or 81 km/h. a — 0.2 corresponds to a maximum 
acceleration of 1.5 m/s per second or 5.4 km/h per second, b — 0.6 corresponds to a maximum deceleration of 
16.2 km/h per second. 

All simulations are done in a 1-lane system of length L with periodic boundary conditions (i.e. the road is 
bent into a ring). Let N be the number of cars on the road. The (global) density is p = . 



3.2 Pictures 



Before analysing the Krauss-model numerically, it is instructive to look at the space-time plots in Fig. 1(c). Space-time 



plots are pictures of the time evolution of the system. In Fig. 1(c) , vehicles drive to the right and time goes down. Each 
row of pixels is a "snapshot" of the state of the road. In principle, one can reconstruct the trajectory of a particular 
car by connecting the corresponding pixels. At the displayed resolution this is however close to impossible and it is 
mostly the larger scale traffic jam structure that one observes. Traffic jams move against the direction of driving. 

i) The laminar state: All cars drive at high speed. The available space is shared evenly among the cars. The 
traffic is very homogeneous. 

ii) The mixed state: The slow cars are all together in one big jam. On the rest of the road, the cars drive at high 
speed. The traffic is very inhomogeneous. 

iii) The jammed state: The density is so high that not a single car can drive fast. As in a), the traffic is very 
homogeneous. 

iv) The single phase at high e: Many small jams are distributed over the whole system. There is neither a larger 
area of free flow, nor a major jam. The traffic is homogeneous. 

Note that "homogeneous" here means "homogeneous on large scales". What this means is that there is a 
spatial measurement length £ above which all density measurements return the same value.[] If a system goes 
from a 2-phase to a 1 -phase model, then even in the regime which technically has only one phase, structure 



formation on small scales is still possible. Fig. 1(c) bottom right is indeed an example for this. Only with 



larger distance from the 2-phase model will these structures go away. 



3.3 Defining a jam 

In order to make progress, one needs to define where a jam starts and where it ends. Our definition of homogeneity 
(see later) will not depend on this, however A jam is a sequence of cars driving with speed less or equal Umax/2. The 
cars between two neighbouring jams are in laminar flow. 

This definition is very simple, but will not always correspond to our natural understanding of the word jam. 
Thus, whether a car is jammed or not according to this definition is just a starting point and not the final answer. 

' More precisely, the fluctuations from one measurement to the next are the same as in case (i). 
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3.4 Initial Condition And Relaxation 

For many parameters of the Krauss model, there is a unique equilibrium state, which the system will attain after a finite 
time freiax, HO matter how it was started. However, deciding when the equilibrium is reached is not trivial (running the 
simulation for t ^ oo clearly is not an option). 

Let Tt be the state of the road at time t. To find the equilibrium value of some property, we use the 

following idea: For small t, E[f{rt)] will depend on the initial condition. With increasing time, E[f{rt)] converges 
towards the equilibrium value. Assume the convergence is from above. Now we need another initial condition that 
approaches the equilibrium value from below. Once these two sequences are close enough together, the equilibrium 
value is found. Unfortunately, it cannot be guaranteed that the value thus obtained really is the equilibrium value. 

We use the following two initial conditions: 

• laminar: The cars are positioned equidistant over the road with speed zero. 

• jammed: all cars are cramped together in a big jam without any gap. Their speed is zero. 

An example of this method is shown in Fig. 0. For /(.) the number of jams was used. Since both initial 
conditions start with t; = 0, the criterion finds one large jam. Vehicles then accelerate, but because of interaction will 
form small jams. For that reason, the laminar start leads to many jams very quickly. From then on, the number of jams 
goes down, because jams coagulate. In contrast, when starting with a large single jam, than that jam remains the only 
one in the system for large times. In Fig. ^ we see that for e = 1.0, the system eventually goes to a state where it has, 
in the average, about 1.8 jams. In contrast, with e — 1.5, the system converges to an average of more than 20 jams. 
Also, the figure shows that the system goes to those long-run states no matter how it starts. 



4 ESTABLISHMENT OF A PHASE DIAGRAM VIA A MEASURE OF IN- 
HOMOGENEITY 

In this section, a criterion is established that distinguishes homogeneous from mixed states. As pointed out before, 
mixed states are characterized by the coexistence of laminar and jammed traffic. As we have also said, if we wait long 
enough, the phases will coagulate, leading to exactly one laminar and one jammed section in the system. When 
approaching the boundaries of the mixed regime, this characterization will become less clear-cut, and it may be 
possible to have more than one jam. Typically, there will be one major jam and many small ones, but for many 
measurement criteria this will cause enough problems to no longer be able to differentiate between the mixed and a 
homogeneous state. Our criterion will only show a gradual decrease in differentiating power. 

As already discussed earlier, it should be noted that some states that are called "homogeneous" in this paper 



may appear inhomogeneous to an observer. An example for this is Fig. 1(c) bottom right. As said before, these states 
are "homogeneous on large scales", which is the important criterion here. Essentially, this means that for system size 
L ^ CO and measurement interval £ ^ oo (but I <C L), all density measurements will eventually return the same 
value. This will not be the case for "mixed" states. 

The criterion is defined as follows: Partition the road into segments of length £ (for simplicity let £ divide 
L without remainder). For each segment the local density pi can be computed as the number of cars in that segment 
divided by i. An interesting value is the variance of the local density: 

L/l 

VarH-— ^(p,(z)-E[p,])2, (4) 

' i—l 

where E[] is the expected value, which in our case is the same as the systemwide density. Note that since the density 
lies within [0, 1], the variance cannot exceed 1/4. 

What this value picks up is how much each individual measurement segment of length £ deviates, in terms 
of its density, from the average density. Assume a system consisting of jammed and laminar traffic. If there is a jam 
in one segment, then the segment's density will be much higher than the average density. Conversely, if there is only 
laminar traffic in a segment, then the segment's density will be much lower than the average density. Var[p£] takes the 
average over the square of these deviations. 

Fig. 1^ shows this value as a function of the global density p and the noise parameter e. Look at it for fixed e, 
say e = 1. One sees that at densities up to p « 0.2, the value of Var[p£] is close to zero, indicating a homogeneous 
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state, which is in this case the laminar state. Similarly, for densities above p k, 0.8, Var[p£] is again close to zero, 
indicating a homogeneous state, which is in this case the jammed state. In between, for 0.2 < p < 0.8, the value of 
Var[pf] is significantly larger than zero, indicating a mixed state. 

Now slowly increase e. We see that the laminar regime ends at smaller and smaller densities, while the 
jammed regime starts at smaller and smaller densities. The latter means that for large e, the jammed phase has many 
relatively small holes, which reduce the density, but do not break the j am. At e w 1.7, the mixed phase completely goes 
away; for larger e, we do not pick up any inhomogeneity at any density. Compare this to the theoretical expectation 
in Fig. 1(b), where for increasing T the two densities eventually merge and thus the different pha ses go away. Note 
that close to the transition the system still looks like it possesses different phases (see Fig. |l(c)| (iv) and locate the 
corresponding e = 1.8 and p = 0.2 in Fig. H). These structures do however exist on small scales only; when averaging 
over larger segments, then all segments contain exactly the same density. A segment lengt h of t = 62.5, as used in the 
figure, is already sufficient in order to not measure any inhomogeneity for the state in Fig. 1(c) ^iv). 

Remember again that e is a model parameter while p is a traffic state. That is, once one has settled for an e, 
the model behavior is fixed, and one has decided if one can encounter a second phase or not. If one can encounter a 
second phase, it will come into existence through changing traffic demand throughout the day - traffic can move from 
the laminar into the mixed and potentially into the jammed state and back. 

As a side remark, let us note that there is also another 1 -phase regime for e 0. Albeit potentially interesting, 
this is outside the scope of this paper. 

The maximum of the density variance is at p « 0.5, e « 1. This can be explained as follows: e = 1 produces 
a sharp separation of the jam (high density) and the laminar structure (low density). With p = 0.5 it turns out that 
these two phases have the same length, thus Var[p^] is maximal. Increasing (decreasing) p will increase (decrease) the 
length of the jam and therefore in both cases decreases Var[/5^]. 

The elliptical shape with its diagonal axis (instead of a vertial axis as in Fig. 1(b)) can be explained as follows: 
Increasing e decreases, via Eq. |^, the acceleration. This leads to lower density in the laminar region. Since Var[p^] 
is maximal when both phases have the same length, and since a lower density in the laminar region makes the jam 
occupy more space, one needs to reduce the density in order to go back to the state where they occupy equal space. 

In summary, one obtains, for the traffic model, a phase diagram as in Fig. 1(b), which is the schematic phase 
diagram for a gas-liquid transition in fluids. Again, the important feature of this phase diagram is that there are three 
states (laminar/gas; mixed/coexistence; jammed/liquid) for low temperatures. For higher temperatures, the coexistence 
range becomes more and more narrow, while the density of the gas phase and the density of the liquid phase approach 
each other. Eventually, they become equal, and the coexistence state dies out. The only important difference is that for 
our traffic model the phase diagram is bent to the left with increasing e. 

There are other criteria which can be used to understand these types of phase transitions. In particular, one 
can look at the gap distribution between jams, and one would expect a fractal structure at the point where the 2-phase 
and the 1-phase model meet, i.e. at p w 0.2 and e « 1.7. This is still under investigation but looks promising. This 
should be followed by an investigation in how far cell-based traffic models ( |J^ models display a similar transition. 



5 PHASE TRANSITIONS IN DETERMINISTIC MODELS 

Only stochastic models can display spontaneous transitions between homogeneous and coexistence (= mixed) states. 
The nature of the transition can however also become clear in deterministic models. We will discuss these similarities 
first for a deterministic car following model and then for deterministic fluid-dynamical models. 

5.1 Car Following Models 

For the model of Eq. (|l]), it has been shown (iS) that the homogeneous solution of the model is linearly stable for all 



densities for V' < a/2, and linearly unstable for certain densities for V' > a/2, where V' is the first derivative of the 
function V{.). The instability sets in for certain intermediate densities, that is, for low and high densities all models 
are stable in the homogeneous (laminar or jammed) state. For intermediate densities, one can select the curve V{g) 
and the parameter a such that the model is either stable or not stable. 

If all parameters including the density are such that the homogeneous solution is not stable, then the system 
rearranges itself into a pattern of stop-and-go traffic, corresponding to the coexistence state. The density of the laminar 
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and the jammed phase in the coexistence state are independent from the average system density, that is, if in that state 
system density goes up, it is reflected in the jammed phase using up a larger fraction of space. 

The type of the instability is similar to the better-known instability of Eq. ([l|). However, once the instability 
is triggered in Eq. (jl]), it will just grow exponentially, and no stable 2-phase solution is found (e.g. ([7^). 

5.2 Fluid-Dynamical Models 

Standard Lighthill-Whitham theory, of the type 

dtp + dxQ{p) = 

with a strictly convex flow-density-curve Q{p), results in a 1-phase model. When Q{p) has linear sections, then in 
that range shock waves are marginally stable, in the sense that disturbances to those shocks are neither amplified nor 
dissipated away. (In a 1-phase model, disturbances are dissipated away and the final state is always homogeneous; 
in a 2-phase/3-state model, there is a density regime where disturbances to shock waves are dissipated away in the 
sense that the shock wave will regain its standard profile again.) When Q{p) has concave sections, then the situation 
becomes more complicated (e.g. (|2^). 

Fluid-dynamical theory, of the type 

dtp + dx{pv) = 

and 

dfV + V dxV 

can, depending on the choice of parameters 
model (|7]). 

As pointed out before, these models are deterministic, so in no situation will these models display stochastic 
transitions. 

6 DISCUSSION 

There is no general agreement if measurements show 1 -phase/1 -state or 2-phase/3-state traffic (or possibly even three 
phases (|22|)). There is some evidence for hysteresis in Germany (]2^, manifesting itself in transitions from high to 
lower flow values at the same density. Hysteresis, which was also found earlier (|2£), is a strong indication for a 2- 
phase model. However, even in Germany, most measurements indicate highly variable traffic at intermediate densities, 
which does not correspond to any clear-cut picture. 

In this context, one should note that a 1-phase model which is close to a 2-phase model would also display 
highly variable traffic at intermediate densities, although it would be homogeneous at large scales as discussed in 
Sec. This variability is however a property of stochastic models only and for that reason it is not well integrated 
into current theory development. A precise investigation of these relations is beyond the scope of this paper. It seems 
however impossible to us to clarify the question if traffic displays several phases or not - and therefore, if breakdown 
probability should be entered into the Highway Capacity Manual or not - without having understood how different 
phases are generated by stochastic models. The present paper fills exactly this gap. 

7 SUMMARY 

This paper shows, via numerical evidence, that a specific stochastic car following model can either display 1 -phase/1 - 
state or 2-phase/3-state traffic, depending on the choice of parameters. With 2-phase parameters, the two phases are 
laminar and jammed, which also corresponds to two of the three states. Those states are homogeneous. The third 
state, at intermediate densities, is a coexistence or mixed state, consisting of sections with jammed and sections with 
laminar traffic. 

The transition to a 1 -phase/1 -state model happens via the densities of the laminar and of the jammed phase 
approaching each other until they become the same. Beyond this point, there is only one homogeneous phase of traffic. 



^-{y{p)-v)+a{p)dxp + v{p)dlv 

including the V {p)-cm:\t, either be a 1-phase/l-state or a 2-phase/3-state 



Jost,Nagel 



8 



Some of these findings can be understood by looking at deterministic models for traffic, either car-following 
or fluid-dynamical. However, the stochastic elements of the transition cannot be explained by deterministic models. 
An important stochastic element is that structure formation and strong variability can also happen in a 1 -phase model 
as long as the parameters are close to the 2-phase model - a deterministic model would converge to a homogeneous 
solution here. 

In our view, it is important to understand this possibility of stochastic models to be in different regimes if 
one considers to enter discussions of traffic breakdown probabilities into the Highway Capacity Manual. If traffic is 
best described by a 1 -phase model, then there is, in our view, no theoretical justification for such probabilities. If, 
however, traffic is best described by a 2-phase model, then the 2-phase model could even give theoretical predictions 
for breakdown probabilities. A discussion of breakdown probabilities in 2-phase models can be found in Ref. (^J). 
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FIGURE 1 : (a) Schematic representation of the gas-Uquid transition in one dimension, (b) State of the gas-fluid 
model as a function of the density and the temperature, (c) Space-time plots for different parameters. Space 
is horizontal; time increases downward; each Une is a snapshot; vehicles move from left to right; fast cars are 
green, slow cars red. L = 600 for all plots. 
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FIGURE 2: Time evolution of the number of jams. All four curves are for 1000 cars. Each curve is an average 
over many realizations, each with a different random seed. 
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FIGURE 3: 3d-plot and isolines of the density variance. The outermost isoUne is Var[p£] = 0.01, the innermost 

Var[pf] = 0.09. L = 1000 and t = 62.5 



